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Abstract 

We propose a description of pairing properties in finite systems within the canonical and micro- 
canonical ensembles. The approach is derived by solving the BCS and self-consistent quasiparticle 



Q . random-phase approximation with the Lipkin-Nogami particle-number projection at zero temper- 

P5 , ature. The obtained eigenvalues are embedded into the canonical and microcanonical ensembles. 



The results obtained are found in quite good agreement with the exact solutions of the doubly- 
folded equidistant multilevel pairing model as well as the experimental data for ^^Fe nucleus. The 



> 

OO 

[^^ , merit of the present approach resides in its simplicity and its application to a wider range of particle 



number, where the exact solution is impracticable. 
PACS numbers: 21.60.-n, 21.60.Jz, 24.60.-k, 24.10.Pa 



The study of thermodynamic properties of finite systems such as atomic nuclei or ultra- 
small metallic grains has been an important topic in nuclear physics. Thermodynamic 
properties are usually described within the grand canonical ensemble (GCE), canonical en- 
semble (CE), or microcanonical ensemble (MCE). The GCE consists of identical systems 
in thermal equilibrium, each of which exchanges its energy and particle number with the 
external reservoir. The systems of the CE exchange only their energies whereas the particle 
number is fixed. The MCE consists of thermally isolated systems with the same energy 
and particle number. Among these ensembles, the GCE is often used in the theoretical 
studies because of convenience in the calculations of thermodynamic quantities. For exam- 
ple, the well-known BCS theory of superconductivity [1] was derived based on the GCE. 
This theory describes very well thermodynamic properties of infinite systems such as metal 
superconductors, where quantal and thermal fluctuations are zero. The latter, however, 
have been shown to be quite significant in finite small systems |2|-l7|. They smoothed out 
the superfluid-normal (SN) phase transition, which is a typical property of infinite systems. 
Most of theoretical approaches to thermal pairing have been derived so far within the GCE 
in finite systems such as atomic nuclei, where no particle-number fluctuations are allowed. 
Therefore, the CE and MCE should be used instead. Moreover, the exact eigenvalues of 



fl&. 



But this task is 



the pairing problem |8| are usually embedded into the CE and MCE 
impracticable for particle numbers A^ > 14 in the case of half-filled doubly-folded multilevel 
model with N = Q {Q is number of single-particle levels). The reason is that the exact 
partition function at finite temperature T 7^ should include all excited states. The finite- 



temperature shell-model Monte Carlo (FTSMMC) 10|, lUj has also been used to evaluate 
the partition function, but it is very time consuming and cannot be applied to heavy nuclei. 
The FTSMMC cannot be applied to the MCE either because it is impossible to include all 
the microstates of the system. It is therefore highly desirable to construct an approach based 
on the CE and MCE, which can offer results in good agreement with the exact solutions for 
any value of the particle number, as well as with experimental data for realistic nuclei. This 
is the goal of the present study. 

We consider the pairing Hamiltonian H = J2j^j^j " ^J2jj'^j^j' > where Nj = 
Y.m^]m^jm IS the particlc-number operator, and Pj = J2rri>o ^]m(^]rh^ Pj = i^jV are the 
pairing operators. The operators a^^ and aj^ are respectively the creation and destruction 
operators of a nucleon moving on the j-th orbital {ajm = (— )"'^'"aj_m)- This Hamiltonian 



a 



has been diagonalized exactly in [8|] to obtain nExact=X]s'C5 ^ ^a^"- -5/2 eigenstates with 
eigenvalues e^^^'^^ and occupation numbers f^, where S = 0, 2, ... Q is the total seniority 
of the system and C^ = m\/[n\{m — n)\] [7]. The exact partition function is constructed 
by embedding the exact eigenvalues into the CE as ^Exact(/3) = Xls '^s'^^P(~/^^i^^'^*) > 
with the degeneracy ds = 2^ and inverse temperature /3 = 1/T. Knowing the partition 
function Z, one calculates the free energy F, entropy S, total energy £, heat capacity C, 
and pairing gap A as F = -rinZ(T), S = -dF/dT, S = F + TS, C = dS/dT, and 
A = [—G{£ — 2 ^ . ejfj + G ^ . ff)]^^"^, where fj are the occupation numbers on the jth or- 
bital obtained by averaging the state-dependent occupation numbers / • within the CE |7| . 
The conventional BCS theory at T 7^ (FTBCS) is derived by using the variational proce- 
dure with respect to the quasiparticle Hamiltonian "H. The latter is obtained by applying the 
Bogoliubov's transformation of the pairing Hamiltonian H from the particle operators, aj^ 
and ajm, to the quasiparticle ones, a]^ and ajm- The variational procedure is carried out 
on the average quantities within the GCE, namely (O) = TrfOe-^^^-^^VTrie-^^^-^^^] 
for any observable O with A being the chemical potential determined as the Lagrangian 
multiplier to preserve the average particle number (A^) = A^. Therefore the conventional 
FTBCS theory is called the GCE-BCS in the present article. 

The CE description of thermal pairing can be undertaken in two directions. The first 
one is to apply the exa ct p article-number projection (PNP) at finite temperature on top 



of the GCE theory [12, ll3|. This approach is rather complicated for realistic nuclei. The 
second direction, which the present study follows, is to embed the solutions of a theoretical 
approximation at zero temperature, which conserves the particle number, into the CE. We 
employ the solutions of the BCS equations with PNP within the Lipkin-Nogami method [iJ] 
(LNBCS) for each total seniority 5* and embed the eigenvalues e:^^^'-'^ obtained at T=0 into 
the CE. The partition function of the so-called CE-LNBCS approach is then given as 

s 
Using it, we can obtain the free energy, entropy, total energy and heat capacity as has 

been discussed above, by replacing ^Exact(/3) with ^lnbcs(/3)- As for the pairing gap, it is 

obtained by averaging the gaps A^^^*-'^, which are the solutions of the LNBCS equations at 

T = 0, namely 

A = Zlnbcs(/3)-^ Y. A^ ^'^56-^^^^''^ . (2) 



Within the BCS theory at T = 0, only the lowest eigenstates can be obtained, e.g the 
ground-state energy for 5=0. The total number of the LNBCS eigenstates embedded into 
the CE is equal to ulnbcs = Xls^^, which is much smaller than UExact- The CE of these 
lowest eigenstates is therefore comparable with the exact one only in the region of low T. 
At higher T, one needs to include not only the ground state but also excited states into 
the CE. This can be resolved by means of the selfconsistent quasiparticle RPA with Lipkin- 



Nogami PNP (LNSCQRPA) [l5J. The LNSCQRPA includes the ground-state and screening 



correlations, which are neglected within the conventional BCS and QRPA. The importance 



of these correlations in finite systems has been demonstrated in 15|. They improve the 
agreement between the energies of ground state and low-lying excited states obtained within 
the LNSCQRPA and the corresponding exact results for the doubly-folded multilevel pairing 



model. The formalism of the LNSCQRPA was presented in details in 15| , so we do not repeat 
it here. The total number of eigenstates obtained within the LNSCQRPA is ulnscqrpa = 
J2s *-"§ X (i7 — S) >nBcs because of the presence of QRPA excited states u, but it is still much 
smaller than UExact- For example, for the half-filled model with Q = N = 10, UExact = 8953, 
whereas ulnbcs = 512 and ulnscqrpa = 256L The thermodynamic quantities are obtained 
within the CE-LNSCQRPA in the same way as that for the CE-LNBCS ([1]), namely from 
the CE-LNSCQRPA partition function Zlnscqrpa(/3) = Es ^5exp[-/3e5^^^^^^^], where 
€g are the eigenvalues obtained by solving the LNSCQRPA equations for each total 

seniority 5*. 

Within the MCE description, we use the eigenvalues ^lnbcs ^^^^ ^^ Q ^^ calculate 
the MCE entropy directly from the Boltzmann's definition S{S) = InW (S) ^heie W{S) is 
the number of accessible states within the energy interval {£,£ + 58) [7|, [l7|. The corre- 
sponding approaches, which embed the LNBCS and LNSCQRPA eigenvalues at T = into 
the MCE, are called the MCE-LNBCS and MCE-LNSCQRPA, respectively. 

The proposed approach is tested within a schematic model and applied to describe ther- 
modynamic quantities for a realistic nucleus. For the schematic model, we employ the 
so-called Richardson (or ladder) model, which consists of Vt doubly-folded equidistant lev- 
els with the single-particle energies chosen as e^ = j — (i7 + l)/2 MeV. The model is 
half-filled, i.e. N = Vt. As for the realistic nucleus, we employ the axially deformed 
^ The first solution of tlie SCQRPA equations corresponds to tire spurious mode and is subtracted from the 
total number of solutions. 
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FIG. 1. (Color online) Pairing gap A, total energy £ and heat capacity C as functions of T 
obtained within the Richardson model for A^ =10 with G=\ MeV. The thin dashed, thick dashed, 
thick sohd, and thin sohd hnes denote the GCE-BCS, CE-LNBCS, CE-LNSCQRPA and exact CE 
results, respectively . 
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FIG. 2. (Color online) Pairing gap A, free energy F, excitation energy £"*, heat capacity C, CE 
and MCE entropies S as functions of T for ^^Fe. In (a), the solid and dash-dotted lines denote the 
neutron and proton gaps, respectively, with the thin and thick lines corresponding to CE-LNBCS 
and CE-LNSCQRPA results, respectively. In (b) - (e), the thin and solid lines denote the CE- 
LNBCS and CE-LNSCQRPA results, whereas the dashed line and the full squares connected by 
solid line stand for the experimental data and the FTSMMC results, respectively. In (f), the MCE 
entropies obtained within the MCE-LNBCS (circles), MCE-LNSCQRPA (triangles) are plotted vs 
experimental data (open squares with error bars) . 



Woods-Saxon single-particle spectra for Fe^^ nucleus with the depth of the central potential 

V = 1/o[l ± k{N - Z)/{N + Z)], where Vq = 51.0 MeV, k = 0.86, the diffuseness a=0.67 



fm, ro=1.27 fm [161], the deformation parameter /?2 = 0.244 lOJ, and the spin-orbit strength 
A = 35.0. 

Shown in Fig. [T]are the pairing gaps A, total energies £ and heat capacities C obtained 
within the GCE-BCS, CE-LNBCS, CE-LNSCQRPA as well as the exact CE for the Richard- 
son model with A^ =10 and G = 1 MeV. The figure clearly shows that the CE-LNSCQRPA 
results (thick solid lines) nearly coincide with the exact CE (thin solid lines) for all thermo- 
dynamic quantities under consideration. A slight difference between the CE-LNSCQRPA 
results and the exact ones at very high T occurs because the CE-LNSCQRPA includes only 
the low-lying excited states in the ensemble. These states are known to be important in low 
and medium temperature regions. At very high T, one needs to include high-lying excited 
states into the ensemble as well. The same feature is seen from the results obtained for other 
systems with different A^ and G. The results obtained within the CE-LNBCS (thick dashed 
lines) are a bit off from the exact ones but, as compared to the predictions by the GCE-BCS 
(thin dashed lines), they still offer much better agreement with the exact solutions. The 
most interesting feature is that neither the pairing gaps obtained within the CE-LNBCS nor 
those obtained within the CE-LNSCQRPA collapse at the critical temperature Tq as pre- 
dicted by the GCE-BCS, but they all monotonously decrease with increasing T, just like the 
exact CE gap. Consequently, the sharp peak in the heat capacity, which is the signature of 
SN phase transition, is smoothed out within these approaches. This feature implies that the 
CE-LNBCS and CE-LNSCQRPA include the effects of quantal and thermal fluctuations, 
which are neglected in the GCE-BCS. Moreover, it also suggests that, if any exact PNP 
method at flnite T could bring the GCE-BCS to the CE-BCS, the results obtained within 
this projected CE-BCS should have the same qualitative behavior as the CE-LNBCS results 
presented in present article with a nonvanishing gap. 

The pairing gap, free energy, excitation energy, heat capacity and entropy in ^^Fe nucleus 
are displayed in Fig. El In order to have a consistent comparison with the recent experimental 
data in |17| as well as the FTSMMC results for these quantities, we employ the same 
configurations proposed in [Uj , namely the CE-LNBCS and CE-LNSCQRPA are calculated 
within the complete pf + OgQ/2 major shell. The results obtained within this shell are then 
combined with those obtained within the independent-particle model (IPM) by using Eq. 



(15) of [ll|, where Z'^^^^ in the present case is the CE-LNBCS or CE-LNSCQRPA partition 
function, whereas Z^^ and Z ^^ are the CE partition functions obtained within the IPM for 
full (from bottom to /iii/2 orbital) and truncated (from /7/2 to /iii/2 orbitals) model spaces, 
respectively (The width of resonance single-particle states is neglected for simplicity). The 
parameter G is adjusted so that the pairing gaps for protons (Z) and neutrons (N) obtained 
within the LNBCS at T = reproduce the experimental values, namely A^ ^1.57 MeV and 
Ajv ^ 1.36 MeV (l7l. It is clear from this figure that the CE-LNSCQRPA results agree quite 
well with the experimental data, which are also deducted from the CE. The small difference 
between the CE-LNSCQRPA results and experimental data at high T is probably due to 
the fact that latter have been extracted from the experimental level density rather than 
being obtained from direct measurements. The level density, however, cannot be measured 
at high excitation energy E* > 10 MeV, but has been instead obtained from the back- 
shifted Fermi gas (BSFG) formula (see e.g. Eq. (2) of [iTJ), assuming a zero pairing gap 
in this region of excitation energy. As we can see now, this approximation is rather crude 
because the pairing gap is always finite even at high E* > 10 MeV (or T > Tc) because of 
thermal fluctuations. This is demonstrated in Fig. [2] (a), which clearly shows that, at T > 
1.5 MeV, which corresponds to E* > 10 MeV, the neutron and proton gaps are still finite. 
This feature puts the use of BSFG formula to extrapolate the measured level density to 
high E* under question. As regards the heat capacity obtained from the experimental data 
plus BSFG [Fig. [2] (e)], it only shows an oscillation rather than a peak or even an S shape 
in the temperature region around Tc, whereas all CE-LNBCS, CE-LNSCQRPA as well as 
FTSMMC results predict a bump in this region. The absence of a bump at T ~ T^ in 
the experimentally extracted heat capacity for ^^Fe also seems to contradict those obtained 
for nuclei in rare-earth region by using the same method, where a clear S shape or even 



a broad bump has been seen 18|. Regarding the difference shown in Fig. [2] (e) between 
the FTSMMC results and ours at high T, it might come from the absence of quadrupole 
interaction in our calculations. It has been well-known that the multipole interactions, in 
particular the quadrupole force, which is important for quadrupole deformation, can have 
significant effects on the nuclear level density at moderate and high temperatures (See Fig. 
8 and Table 1 in Ref. [l9|). Fluctuations of other degrees of freedom such as nuclear shapes 
due to rotation may also play an important role. Most probably these effects will show up at 
high temperatures and angular momenta, which are considered neither in the experimental 



data for Fe that we use nor within the present framework of our modeL 

A genuine thermodynamic observable is the MCE entropy at low E* because it is calcu- 
lated directly from the measured level density (See e.g. Eq. (3) of [17]). Shown in Fig. [2] 
(f) are the MCE entropies obtained within the MCE-LNBCS and MCE-LNSCQRPA using 
the Boltzmann's definition with 6S = 1 MeV versus the experimental data. It is seen that 
the MCE-LNSCQRPA entropy not only offers the best fit to the experimental data but also 
predicts the results up to higher E* > 10 MeV. 

In summary, we have proposed two approximations, which embed the solutions of the 
BCS and SCQRPA with Lipkin-Nogami PNP at T = into the CE and MCE. The proposed 
approaches are tested within the Richardson model and applied to describe the recent exper- 
imentally extracted thermodynamic quantities for ^^Fe nucleus. The analysis of numerical 
calculations for the pairing gap, free energy, total energy, heat capacity and entropy shows 
that the CE-LNSCQRPA predictions are in quite good agreements with the exact results 
(whenever the latter are available), as well as the recent experimental data. Moreover, this 
is a microscopic approach that can describe simultaneously and selfconistently all the ex- 
perimentally extracted quantities, namely the free energy, total energy, heat capacity, and 
entropy within both CE and MCE treatments. To our knowledge, this is the first time that 



the recent experimental MCE entropy 17| has been successfully described by a consistent 
microscopic theory. It also shows that the SN phase transition predicted by the conven- 
tional GCE-BCS theory is smoothed out even within the CE-LNBCS calculations due to 
the effects of quantal and thermal fluctuations, leading a nonvanishing pairing gap. The 
results of the present study put under question the use of BSFG formula to extrapolate the 
measured level density to high excitation energy. The merit of the present approach resides 
in its simplicity and its feasibility in the application to larger finite systems, where the exact 
matrix diagonalization and/or solving the Richardson equation are impracticable to find all 
eigenvalues, whereas the FTSMMC method is time consuming. 

The authors thank Peter Schuck (Orsay) for fruitful discussions. NQH is a Nishina Memo- 
rial Fellow at RIKEN. The numerical calculations were carried out using the FORTRAN 
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